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O MÉTODO DE AJUSTAMENTO COMPLETO 
PARA O CÁLCULO DE BARRAGENS ABÓBADA 


RESUMO 


No metodo de ajustamento completo considera-se à bar- 
ragem dividida em dois sistemas de elementos — arcos e con- 
solas — que são formados por aduelas, Nas aduelas de eru- 
zamento de arcos e consolas são consideradas todas as for- 
cas, internas e externas, e estabelece-se o conjunto completo 
das equações de equilibrio. São consideradas, também, 
as três componentes do deslocamento linear e as seis rotações 
sofridas pelas normais às faces de cada adueta, Mostra-se 
que, por razões de continuidade e equilibrio das tensões 
de corte, as torções segundo os eixos horizontal e vertical 
não são independentes e que os deslocamentos horizontais 
tangenciais e os verticais devem ser considerados simul- 
tanecamente 

As forças internas são obtidas pela resolução de um sis- 
tema de equações de compatibilidade que exprimem a igual” 
dade dos deslocamentos em todas as aduelas de cruzamento 
de arcos e consolas. Ambos os sistemas de elementos são ana- 
lisados para forças triangulares unitárias, para as quais 
se determinam os respectivos coeficientes de influência, Uti- 
liza-se uma formulação matricial para v cálculo dos arcos 
e uma fórmula geral para as consolas, O método pode ser 
usado em computadores de capacidade média, 

Os resultados obtidos pela aplicação do método a duas 
barragens tipo, foram comparados com: os obtidos por 
mésodos britânicos recentemente desenvolvidos, 

Apresenta-se também comparações com resultades de 
ensaios em modelo, A concordância é excelente, 


o. NOTAÇÃO 


Usam-se os seguintes símbolos neste trabalho : 


Ego: St 
D,, 


== matrizes que exprimem os 
deslocamentos devidos às 
cargas actuantes na aduela ; 
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The method considers the arch dam divided in tao sys- 
tems of elements — arches and cantilevers — which are form- 
ed by voussorrs, In the voussoirs of crossing of arches and 


cantilevers all the internal forces and external loads are 
considered and the complete set of equatrons of equibibrium 


is established, The three components of linear displacementsa 
and the six rotations of the normais to the faces are con- 
sidered, For reasons of continuity and of equilibrium of shear 
stresses it is shown that the torsions of horizontal and ver- 
tical axes are not independent and that the vertical and 
tangential displacementa must be taken simultaneously, 

The internal forces are obtained by solving a system 
of equations of compatibility, expressing equal displacement 
at the crossing of arches and cantilevers, Both systems 
of elements are analysed for triangular wnit loads, for 
which the coefficients of influence of the various displace- 
ments are determined, A matrix formulation for the arches 
and a general formula for the cantilevers are used, The 
method can be used in medium size computers, 

The results obtained by the application of the method 
to tuo test dams are compared with those of recently devel- 
oped British methoda. Also comparisons with results of model 
testa are presented, The agreement obtained is excellent. 


B,; == matriz de transformação dos 
esforços unitários na secção 
j de um arco em esforços 
na secção i; 

d = distância do centro de gravi- 

dade ao paramento de mon- 

tante na secção j da consola; 
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k 


N,N,, No, N, 
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= módulo de elasticidade 
(Young) e módulo de elas- 
ticidade transversal; 


= matriz que exprime as for- 
ças devidas às cargas uni- 
tárias distribuídas ; 


= espessura na secção j; 


= funções definidas pela Eq 45 
usada na fórmula para os 
deslocamentos das consolas ; 


== momentos de inércia; 


== tangente do ângulo entre a 
linha dos centros das con- 
solas e a direcção vertical; 


= relação de distorção devida 
ao corte (normalmente 
K=1,5): 


== coeficientes de influência 
das cargas unitárias nos 
deslocamentos dos arcos ou 
das consolas; 


= coeficientes que represen- 
tam a grandeza na aduela 
da carga triangular unitária ; 


= momentos flectores; 


== momento flector na base da 
consola devido à subpressão; 


== momentos de torção; 


= momentos devidos ao peso 
e à componente vertical de 
pressão hidrostática ; 


= forças de flexão derivadas 
do momento de torção M,; 

= força de flexão definida pela 
Eg. 10; 


== forças internas normais; 


= impulsão na base da con- 


bh 


sola devida à subpressão; 


== forças de corte ou esforços 
transversos ; 


Q,0,,0, 


== esforços transversos radiais 
ou forças de corte internas 
radiais ; 


= peso da barragem mais com- 
ponente vertical da pressão 
hidrostática ; 


= componente horizontal da 
pressão hidrostática ; 


Ja = (g—g,) =forças fictícias radiais para 


R 

Ryo Ri 

5, 

teot, 

U, 

U,V,W 
Vas Var Wa 
Uç, Veo Wa 
Es É 

Y; 

FA 


arcos e consolas; 


= raio horizontal de curva- 
tura da superfície média; 


= raio horizontal no para- 
mento de montante ; 


== área da secção j das conso- 
las; 


= derivadas das forças de 
corte ou forças internas 
normais ; 


= carga triangular unitária na 
secção j; 


== deslocamentos tangenciais, 
verticais e radiais; 


== deslocamentos dos arcos; 
== deslocamentos das consolas; 


== coordenadas rectangulares; 


== distância vertical desde o 
coroamento até à secção j; 


== distância da face de mon- 


bo 


tante à superfície de refe- 
rência ; 


GR e * o . . 
Los Pos Tor or p= Características elásticas da 


BB, 


fundação definidas pela 
Eq. 38; 


== deformações angulares de- 
vidas ao momento flector 
entre as faces verticais dos 
arcos e as faces horizontais 
das consolas ; 


== distorções entre as faces 
das aduelas; 
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dx 9, = distorções radiais por corte 
nas faces horizontais e ver- 
ticais da aduela; 

tao Es == deformações lineares entre 
duas faces opostas de uma 
aduela; 

7 == distorção devida ao esforço 
transverso definido pela 
Eq. 8; 

Pe do = rotações de eixo vertical 
dos arcos e das consolas; 

Ú = coeficiente Poisson ; 

Cais Cas Pê = tensões tangenciais no para- 


mento de montante, no 
paramento de jusante e no 
centro da secção da aduela ; 


Yam = rotação do eixo nn em torno 
do eixo mm; 

6, 4, == ângulo horizontal entre o 
eixo da aduela e o eixo da 
barragem ; 

O, By Dy == torções das faces da aduela. 
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Número de barragens acma de 500 pes 


1. INTRODUÇÃO 


Nas condições económicas actualmente exis- 
tentes é cada vez maior o interesse pela constru- 
ção de grandes barragens para rega, aproveita- 
mentos hidroeléctricos, abastecimentos de água 
e controle de cheias. São muitas as regiões pobres 
do mundo que podem aspirar a desenvolver-se 
se os seus rios forem regulados e controladas 
as suas cheias. 

Ainda que nos países mais desenvolvidos 
se estejam actualmente construindo menos bar- 
ragens, em grande parte porque os melhores 
sítios já foram usados e a energia hidroeléctrica 
está explorada quase completamente, verifica-se 
um incremento na construção de barragens 
em todo o mundo. A estatística das barragens 
muito altas existentes (acima de 150 metros 
de altura), que se apresentam na Fig. 1, demons- 
tra que as barragens de betão, e muito especial- 
mente dos tipos em abóbada e em casca, têm 
sido as usualmente escolhidas para vales pro- 
fundos e para muito grandes empreendimentos *. 
O problema estático destas barragens é muito 
complexo e a determinação teórica do campo 
de tensões nestas estruturas e nas suas funda- 
ções tem sido sempre considerado pelos enge- 
nheiros civis como um dos problemas mais inte- 
ressantes da teoria da elasticidade e da teoria 
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Fig. 1 — Número e altura de Barragens acima de 500 pés (150 m) por décadas 
(segundo o Registo Mundial de Barragens ICOLD — 1965) 
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das cascas. Nas últimas décadas (especialmente 
nos anos cinquenta) este facto deu lugar ao uso 
muito generalizado dos modelos para ajudar 
os engenheiros a conhecer correctamente o campo 
das tensões nas estruturas que projectam *. Com 
o desenvolvimento dos computadores, não só os 
métodos matemáticos existentes, que exigem 
um trabalho numérico muito longo, se tornaram 
mais divulgados e aperfeiçoados, como também 
em muitos centros de investigação de todo 
o mundo foram imaginados e postos em prática 
novos métodos para o cálculo de barragens arco 
com bases completamente novas. Estudos, pre- 
viamente programados e organizados em vários 
países, e em especial nas universidades britã- 
nicas”, deram como resultado que se obteve, 
na engenharia civil, algo de bastante importante 
e único — que é poder dispor de métodos exactos 
e completos para a determinação de tensões nas 
barragens abóbada e nas suas fundações. 

Estes métodos dividem-se em três categorias: 
(1) métodos de arcos-consolas (dos quais o mais 
desenvolvido foi o «trial.load» *); (2) métodos 
baseados na teoria das cascas (o mais conhecido 
é o método de Tolke); (3) métodos numéricos 
baseados na teoria da elasticidade tridimensional 
(o mais notável destes últimos foi o método 
da relaxação dinâmica, aplicado ao cálculo da bar- 
ragem de Donkan *). Pode dizer-se que em 1968; 
como consequência do aperfeiçoamento de mui- 
tos destes métodos, se podem obter valores coinci- 
dentes quer para as tensões quer para os deslo- 
camentos. A maior parte destes métodos pode 
também ser usada para barragens de dupla curva- 
tura de qualquer forma, apoiadas em fundações 
deformáveis. 

Este trabalho indica as bases de um método 
que realiza a síntese entre o primeiro e o segundo 
grupo atrás referidos, isto é, recorre ao artifício 
de analisar «arcos» e «consolas», mas como por- 
ções da casca que efectivamente é a forma indi- 
visível da barragem. O método inspira-se nas 
equacões da Teoria das Cascas que, mesmo 
quando se analisam separadamente «arcos» e «con- 
solas», respeita com grande aproximação. São 
absolutamente satisfeitas as equações de equilí- 
brio considerando todas as tensões integradas 
em esforços (esforços normais, esforços transver- 
sos, momentos flectores e de torção) como simpli- 
ficações aceites na Resistência dos Materiais 
e recorre-se a equações de compatibilidade, algu- 
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mas sob a forma de ajustamentos entre «arcos» 
e «consolas». Chamou-se-lhe «Método de Ajus- 
tamento Completo» em réplica ao «trial load» 
e a outros métodos de ajustamento que o não 
são completamente, mas é também um método 
de resolução das equações diferenciais da Teoria 
das Cascas, embora não por integração directa 
de todas as equações, como tem sido tentado 
por vários autores. 

Apresentam-se ao longo do trabalho compa- 
rações com resultados de tensões e deslocamen- 
tos obtidos por outros métodos (especialmente 
os de terceiro grupo) ultimamente apresentados 
no Simpósio de Barragens-Arco, realizado em Lon- 
dres”. Também se apresentam resultados da apli- 
cação do método de ajustamento completo a algu- 
mas barragens recentemente projectadas, ao mesmo 
tempo que os valores obtidos em ensaios de mode- 
los dessas mesmas barragens. Como se verá, 
os resultados dessas comparações apresentam 
concordância notável. 

Ainda que este método tenha sido desenvolvido 
para barragens abóbada e em casca, ele pode 
ser usado para o estudo de qualquer corpo 
contínuo delgado analisando um número restrito 
de elementos discretos, e pode ser aplicado ao 
estudo de cascas de forma complexa com espes- 
suras variáveis. O método tem a vantagem de 
ser fácil de aplicar e de permitir o uso de 
computadores de capacidade média. 


2— O MÉTODO DE AJUSTAMENTO 
COMPLETO 


O método baseia-se antes de tudo na Teoria 
das Cascas e, depois, no conceito muito comum- 
mente usado da divisão da barragem em dois 
sistemas de elementos lineares livres: um sistema 
de elementos horizontais (os arcos) e outro 
sistema de elementos ortogonais aos primeiros 
(as consolas). O ajustamento completo dos 
deslocamentos dos elementos de um sistema com 
os deslocamentos dos outros elementos do outro 
sistema, permitirá, em complemento a outras 
considerações impostas pela Teoria das Cascas, 
uma determinação correcta das tensões principais 
na barragem dentro das aproximações habituais 
na engenharia obtidas através da Resistência 
de Materiais. 

Na análise separada dos arcos e das consolas 
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entram em jogo todos os esforços actuantes na 
casca, de que aqueles elementos, arcos ou consolas, 
são simples fatias que justapostas a integram. 
No cálculo dos deslocamentos a ajustar (nem 
todos são calculados, por ser desnecessário como 
se verá) considera-se a contribuição de todos os 
esforços internos da casca. Os esforços conside- 
rados nas secções e nas faces justapostas dos 
«arcos» e das «consolas» são todos os esforços 
internos da casca. 

Nas aplicações da Teoria das Cascas é frequente 
admitir-se, como simplificação, que os segmentos 
de recta normais à superfície média da casca se 
conservam normais e de comprimento indeformado 
(hipótese de Kirchhoff). Com as espessuras 
atingidas nas barragens não é mais admissível 
tal simplificação, e o presente método não a 
admite senão na conservação do comprimento 
indeformável, admitindo que a pressão hidros- 
tática não açtua no paramento de montante, mas 
sim na superfície média. 

Deixou-se para ulterior desenvolvimento a 
consideração do efeito de Poisson, que convirá 
considerar no apuramento do método. 


BARRAGE rs 


Fig. 2 — Elementos arco e consola. Deslocamentos 
e rotações de uma aduela 
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Os arcos são elementos horizontais apoiados 
numa fundação deformável. Eles são definidos 
(Fig. 2) como uma porção de casca compreendida 
entre dois planos horizontais separados por uma 
distância muito pequena (e unitária) na direcção 
vertical. Esta definição, ainda que não seja 
perfeitamente adaptável a barragens arco ou 
abóbada de dupla curvatura, simplifica o problema 
geométrico uma vez que as formas destas 
barragens são definidas por linhas de nível. As 
consolas são elementos ortogonais aos arcos, 
também apoiados numa fundação deformável. 
Quando a curvatura da barragem é diferente a 
cada nível, as consolas são elementos torcidos 
que podem ser gerados por segmentos horizontais 
sempre perpendiculares às linhas de nível da 
face de montante e seguindo nessa face duas 
linhas directrizes ortogonais aos arcos distancia- 
dos duma pequena distância (e unitária) medida 
no coroamento da barragem. 

Para barragens abóboda de dupla curvatura, os 
arcos e as consolas poderiam ser tomados como 
elementos do mesmo tipo, isto é, como elementos 
torcidos normais ao paramento de montante. Os 
arcos seriam então um semi-elemento, indo desde 
cada um dos encontros até ao fecho e seriam 
definidos por uma porção da casca compreendida 
entre as duas superfícies geradas por segmentos 
sempre normais ao paramento de montante e 
seguindo duas linhas de nível afastadas na 
vertical duma distância muito pequena (e unitária). 

Os arcos ou as consolas podem ser considerados 
como uma sucessão de aduelas, dispostas em 
fiadas horizontais ou em fiadas verticais, corres- 
pondendo cada aduela à intersecção de um arco 
com uma consola, quando fossem consideradas 
sucessões contínuas de elementos, arcos e consolas. 
Na prática não é necessário considerar uma 
sucessão contínua de elementos mas apenas 
alguns arcos e algumas consolas. 

O «Método trial load», tal como foi desenvol- 
vido pelo Bureau of Reclamation, toma em 
consideração as componentes radiais e tangenciais 
horizontais do deslocamento linear de pontos da 
superfície média da barragem e rotações de eixo 
vertical do segmento montante-jusante. Nesse 
método desprezam-se as componentes verticais 
do deslocamento, devendo também ser notado 
que as rotações de eixo horizontal-tangencial 
não precisam de ser calculadas. Em contrapartida, 
no método agora apresentado, os movimentos 
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verticais são considerados completemente, assim 
como as rotações de eixo horizontal radial. 
Todavia, uma análise mais completa das defor- 
mações e dos deslocamentos revela que, dada a 
conexão entre umas e outras e a sua dependência 
em relação às tensões, o problema pode ser 
simplificado e o número de incógnitas em cada 
ponto de intersecção dos arcos com as consolas 
(tanto deslocamentos como esforços internos) 
pode ser reduzido. 


3 — ANÁLISE DOS DESCOLAMENTOS E 
DEFORMAÇÕES 


Na Fig. 2 indica-se o sistema de eixos asso- 

ciados com uma aduela. Quando a barragem se 
deforma por acção das cargas exteriores, o centro 
de gravidade de cada aduela sofre um desloca- 
mento linear de três componentes U, Ve W e, 
os três eixos perpendiculares normais aos três 
pares de faces da aduela, sofrem rotações que são 
expressas por seis componentes : 
Vox € Y,y rotações do eixo dos zz em torno do 
eixo dos xx e do eixo dos yy, Y,, € 'V,, rotações 
do eixo dos yy em torno do eixo dos xx e do eixo 
dos zz, v,, € Y,y rotações do eixo dos xx em 
torno do eixo dos yy e do eixo dos zz. 

Em consequência dos esforços do corte, os três 
eixos não permanecem perpendiculares depois da 
deformação e, portanto, as rotações de dois eixos 
perpendiculares em torno de um terceiro eixo 
não são iguais, por exemplo, am + Yax (Fig. 2). 

Das nove componentes dos deslocamentos, 
cinco são independentes, concretamente : os três 
deslocamentos,U, Ve W e as rotações do eixo 
dos zz Y,y e Y,y respectivamente em torno do 
eixo dos xx e do eixo dos yy. As outras rotações 
são derivadas das três componentes dos deslo- 
camentos lineares : 


CC 0W 3 OU 
yx dy ' YZ dy á 
1 90W 4, OV 
Txy dx 1 Txz dx (1) 


As rotações do eixo dos zz 


d 
de Sã 


dy dy 


AP 4 
2x € Yyy diferem 


por causa das distorções 2, e 9, 
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devidas ao corte radial nas faces horizontais e 
verticais da aduela 


Se as distorções º?, e õ, fossem iguais a zero, 
as seis rotações dos eixos deveriam ser iguais 
duas a duas. Neste caso, seria suficiente tomar 
as três componentes dos deslocamentos U, V e 
W para descrever todas as deformações de casca, 
e corresponderia a aceitar a hipótese corrente na 
Teoria das Cascas de que os segmentos normais 
a superfície média permanecem normais a essa 
superfície depois da deformação. Neste caso, em 
que os arcos foram definidos como planos e hori- 
zontais, corresponderia a admitir que os ângulos 
entre os segmentos horizontais normais aos arcos 
a montante e a superfície média permanecem 
constantes. Ver-se-á, no cap. 11, que os erros que 
resultarão por vezes do facto de se desprezar as 
distorções 2, e à, são pequenos na operação do 
ajustamento especialmente nas barragens finas. 
Deve acentuar-se que, quando se calculam os 
arcos e as consolas um a um, estas distorções 
são consideradas no cálculo dos respectivos 
deslocamentos, porque já não é desprezável esta 
contribuição, embora as extensões e as flexões 
permaneçam como principais contribuitnes dos 
deslocamentos. 

Aceitam-se as duas hipóteses básicas da Resis- 
tência dos Materiais: (1) que as faces horizontais 
e verticais da aduela permaneçem planas depois 
da deformação (hipótese de Navier), o que quer 
dizer que as tensões normais e as tensões tangen- 
ciais de montante a jusante, têm uma distribuição 
linear e (2) que são desprezáveis as deformações 
da espessura, correspondentes às tensões normais 
radiais (que são iguais à pressão hidrostática no 
paramento de montante), o que equivale a admitir 
que as cargas são concentradas na superfície média. 

Deste modo, as deformações das aduelas que 
se tomam em consideração no Método, são as 
seguintes, ou as suas combinações, todas derivá- 
veis das componentes dos deslocamentos atrás 
referidos : 


1 — duas extensões ou encurtamentos da 


aduela, E — uma vertical, gs , que 
y 
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é a extensão, na vertical, da linha média 


') 
das consolas, e outra horizontal, *, = — À 
x 
que é a extensão na linha média dos 


arcos ; 


2 — duas flexões É —que são as variações 
angulares entre as faces verticais dos 


OW ey dE 
arcos, É, = e as faces horizontais 
x 
04 
das consolas, É, = —E 
dy 


3 — duas torções w — uma de eixo vertical, 1, 
das consolas, e a outra de eixo horizontal 
tangencial dos arcos, m»,; 


4 — três escorregamentos, as distorções entre 
as faces verticais na direcção radial, 0, 
(corte radial nos arcos), entre as faces 
horizontais na direcção radial 2,, (corte 
radial nas consolas); e distorção angular na 
superfície média y, (distorção tangencial). 
Assim a deformação da linha média dos arcos 
pode ser expressa por: 


— (3) 


em que É, será a deformação atrás definida 
e à, introduzirá o efeito da distorção angular 
devida ao corte radial nas faces verticais. 

Para as consolas, a deformação da linha média 
será expressa também por: 


' 2 
É. + die Vga ow (4) 
Do dy dy dy” 


sendo, analogamente, É, a deformação já atrás 
definida e 2, o efeito da distorção angular devida 
ao corte radial nas faces horizontais. 

Sendo W o deslocamento radial do centro 
de gravidade da aduela, a torção das faces hori- 
zontais , (entre segmentos radiais das faces hori- 
zontais, secções das consolas) é dada por: 


4 2 A 
O Pay e od W/ d e” 


ÉA, + (5) 
dy dxdy dy 


(O = 


y 
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e a torção », das faces verticais (secções dos arcos) 
é dada por: 


W dd, 
o dy Ox 0x 


(6) 


Se aqui forem desprezadas as derivadas das 
distorções angulares em relação às derivadas das 
rotações (Eq. 2), e se se tomar em conta que, 
por razões de continuidade de uma superfície 
deformada, as duas derivadas cruzadas dos des- 
locamentos radiais são iguais, então pode-se 
aceitar que: 


Isto significa que as torções , e V, podem ser 
consideradas iguais tanto para arcos como para 
consolas, e que apenas é necessário considerar 
uma delas, preferivelmente a que se refere 
às rotações em torno do eixo vertical, w,, porque 
é mais fácil de determinar. 

Finalmente, a distorção y, na superfície média 
também se pode escrever em função das compo- 
nentes dos deslocamentos : 


Va dá Vue + Vice So Ea E E EST (8) 


4 — EQUAÇÕES DE EQUILÍBRIO 


No presente método considera-se que as forças 
actuando nas faces radiais de cada aduela repre- 
sentam os integrais das tensões instaladas naque- 
las faces e que estão relacionadas com deforma- 
ções, que por sua vez integram-se, dando os des- 
locamentos, os quais devem ser os mesmos quando 
as deformações são integradas ao longo dos arcos 
ou ao longo das consolas. Deve ainda ter-se 
em atenção que cada um dos sistemas, de arcos 
ou consolas, suporta não só as suas forças inter- 
nas mas também as forças actuando nas suas 
faces limites e que, consequentemente, arcos 
ou consolas representam completamente a estru- 
tura, i.e. cada sistema resiste, por si só, total- 
mente, às forças externas. As forças internas dos 
arcos actuam nas faces limites das consolas 
e vice-versa. 
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Fig. 3 — Forças internas e exteriores actuando numa aduela 


Na Fig. 3 indicam-se as forças internas e as car- 
gas exteriores actuando numa aduela comum 
a um arco e a uma consola. As forças externas 
são o peso mais a componente vertical da pres- 
são hidrostática, representadas por p, o momento 
devido a estas cargas verticais, m, e a compo- 
nente horizontal da pressão hidrostática, q. 

Os métodos simplificados usados no passado 
desprezavam todos eles algumas das forças inter- 
nas indicadas na Fig. 3. Quanto maior o número 
de forças internas considerado, tanto mais exacto 
o método poderia ser. Por exemplo, na «fórmula 
dos tubos» apenas se considerava a força 
interna N,; no método de arcos independentes 
Ne M,e Q,, no método «trial load» de ajus- 
tamento radial, consideravam-se as forças internas 
anteriores e ainda M, e Q; (N, só se conside- 
rava no cálculo das consolas independentes) 
enquanto que no «trial load» completo tomam-se 
em consideração todas as forças atrás referidas 
mais as forças internas M,, e N,,- No presente 
método não só se consideram todas as forças 
internas indicadas na Fig. 3, mas também se con- 
sideram todas as condições de equilíbrio das 
forças actuantes nas aduelas pertencentes simul- 
tâneamente a um arco e a uma consola. 
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As equações de equilíbrio são as seguintes: 


a) radial 


b) vertical 


oN, oN, 
dy dx 


c) tangencial 


9 N 9N Q 
ls E) o RE 2.) Ea den == () (9) 
dx dy R 
d) rotação de eixo horizontal tangencial 


oM, My 
dy d x 


—Q,;—-N,i-m=0 


e) rotação de eixo vertical 
IM, 4 dM ” 
dx dy 


—-Q,+Nyi=o0 
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Aspecto de 2 sectores de ante-distribuidor de uma das 
3 Turbinas Kaplan de 100.000CV destinadas à Central de 
Carrapatelo, da C.P.E., SARL. (Projecto KMW - NOHAB). 


PONTES ROLANTES, GUINDASTES E 
APAR. DE ELEVAÇÃO ESPECIAIS Projecto e fabrico 


TURBINAS HIDRÁULICAS ——— Fabrico segundo licença de A.C. M. de Vevey, S. 4. 
TURBINAS A VAPOR ——— Fabrico segundo licença de Brown Boveri, Cie. 
CALDEIRAS A VAPOR ————— Projecto e fabrico segundo licença de Foster 


Wheeler, Co. 
EQUIPAMENTOS E INSTALAÇÕES 


INDUSTRIAIS 


CONSTRUÇÕES METALOMECANICAS MAGQUE :... 


ALVERÇA DO RIBATEJO- PORTUGAL 
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SE Clio Aa A A À 


v 


> 


SIDERURGIA NACIONAL 
Caves dos trens de laminagem Tl e T2 
Tubagens 


ALTERNADORES (de potência superior a 2000 kVA), COMPENSADORES 
SÍNCRONOS, MOTORES DE GRANDE POTÉÊNCIA . (fabricação parcial). 


EQUIPAMENTOS PARA INDÚSTRIAS QUÍMICAS E PETROQUÍMICAS 


Tubagens, tremonhas, reactores, permutadores, cisternas, etc. 


EQUIPAMENTOS PARA INDÚSTRIAS PETROLÍFERAS 


Torres, balões, permutadores, filtros, aerorefrigerantes, tubagens, etc. 


EQUIPAMENTOS PARA INDUSTRIAS CIMENTEIRAS 

Fornos, moinhos, separadores, arrefecedores, etc. 

RESERVATÓRIOS 

De tecto fixo e tecto flutuante, de alta e baixa pressão. Para gases liquefeitos. Gasómetros 
EQUIPAMENTOS DE MANUTENÇÃO E TRANSPORTE 

Elevadores de alcatruzes, elevadores de palhetas, transportadores de tela metálica e de 
borracha, contentores, ete, 


ESTRUTURAS METALICAS 
Pontes, cimbres, coberturas, pisos, edifícios industriais, etc, 


MONTAGENS PESADAS 


SEPS SOCIEDADE DE CONSTRUÇÕES ELECTRO-MECÂNICAS, S.A.R.L. 


APARTADO 8 + LEÇA do BAILIO+TELEX SEPSA P 2616 
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f) rotação radial 


As equações 9 devem ser satisfeitas em toda a 
barragem. A primeira equação (a) toma em con- 
sideração as forças internas nos arcos e nas 
consolas e, portanto, deve ser verificada em 
ambos os sistemas de elementos. As equações 9 
(c) e (e) serão satisfeitas no cálculo dos arcos e 
as equações 9 (b) e (d) apenas no das consolas. 
A eg. 9 (f) é idênticamente satisfeita por ser 
consequência da igualdade das tensões tangen- 
ciais 7 ,y € 7 yx COMO Se prova. 

Nas equações relativas aos arcos, as forças 
internas nas faces superiores e inferiores podem 
ser substituídas por forças externas estáticamente 
equivalentes, definindo três grupós de forças 
fictícias variáveis q,, t, e m”,. ou seja 


forças radiais 


forças tangenciais 


ON ou 
dy 


== t (10) 


x 


forças de flexão de eixo vertical 


Com estas três variáveis, as três forças inter- 
nas nos arcos N,, OQ. e M, podem ser determi- 
nadas, uma vez que as equações de equilíbrio 
relativas aos arcos são: 


Ne, O, 
R dx 


=. =0 
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oN, Q, 

— — t,=0 Za 
oi R Eee (11) 
oM, 

— OQ. +m',=o0 
dx 


Do que se acaba de dizer, torna-se então claro 
que, neste método, um arco deve suportar não 
só toda a carga exterior da porção de casca que 
ele representa, mas também forças interiores cor- 
respondentes às forças internas nas faces hori- 
zontais que o limitam. 

Também, para as consolas, podem distinguir- 
-se três forças fictícias variáveis: 


forças radiais (as cargas exteriores menos as 
variações das forças internas no arco) 


Ne 00, 
MC 3% =]—J — 
forças verticais 
o N 
— E =t, (12) 
dJ 4 


forças de flexão de eixo horizontal 


d) M sy 


dx 


Das variáveis q. , t, e m, podem determinar- 
-se as forças internas nas consolas, N std e Mo, 
a partir das equações de equilíbrio corresponden- 
tes as consolas, as quais são; 


dQ, | 
E 'q=0 
dy 
N, de 
—p+t,=0 (13) 
y 
dy 
dM 
X -Q,-—-Noi-m+m=0 
0y 
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Torna-se aqui claro, também, que uma consola 
deve suportar não só todas forças externas da 
porção de casca que representa, mas as 
forças internas actuando nas suas faces laterais. 
Pode então concluir-se que o sistema de arcos 
ou de consolas deve, cada um deles, satisfazer 
as equações de equilíbrio no espaço de toda a 
casca, sem desprezar nenhuma das forças internas. 

Conclui-se também que, a partir do conheci- 
mento das forças exteriores, reais e fictícias, é 
possível, analisando cada arco e cada consola 
isoladamente, determinar os esforços e os deslo- 
camentos das secções. 


5 — EFEITOS DOS MOMENTOS DE TORÇÃO 
E ESFORÇOS DE CORTE 


Nos capítulos anteriores escreveram-se as 
equações de equilíbrio e as relações geométricas 
entre deformações e deslocamentos, no que respeita 
às cascas. Falta porém referir as relações entre 
os esforços (ou as tensões) e os deslocamentos 
(ou as deformações), isto é, as relações de 
elasticidade. As equações de equilíbrio agrupa- 
ram-se em duas séries referidas respectivamente 
aos arcos e às consolas, com vista a permitir a 
análise separada desses elementos. A partir dos 
esforços agrupados naquelas séries é efectivamente 
possível determinar os deslocamentos dos arcos 
e das consolas por recorrência às fórmulas bem 
conhecidas da Resistência dos Materiais, relativas 
a estruturas bi-dimensionais. Há porém deforma- 
ções da casca—a distorção tangencial 7, e as 
torções » — que não se podem analisar senão 
visualizando a casca na sua dimensão tridimen- 
sional, deformações essas que têm com os esforços 
relações fundamentais a respeitar. Este capítulo 
dedica-se a estas deformações e às suas relações 
eom os esforços. 

As torções mw, e m, São proporcionais aos 
momentos de torção M,, e M,, que actuam nas 
faces radiais da aduela (Fig. 3). Derivando estes 
momentos de torção obtêm-se forças de flexão 
de eixos horizontais e verticais com as seguintes 
expressões : 


nn, = mn, = —— (14) 


As primeiras forças de flexão actuam nos arcos 
e as segundas nas consolas. Isto é, os momentos 
de torção originam torções nos seus próprios 
eixos e flexões nos elementos ortogonais a esses 
eixos. Estas flexões são consideradas na análise 
separada dos arcos e das consolas. 

No que diz respeito aos deslocamentos tangen- 
ciais e verticais (U e V) da aduela, eles não 
podem ser considerados como independentes, 
como tem sido admitido noutros métodos”. De 
facto, ambos os movimentos estão relacionados 
com as tensões tangenciais médias nas faces 
horizontais e verticais das aduelas 7, e Ty €, 
portanto, devem ser considerados simultânea- 
mente. A distorção por corte tangencial (defor- 
mação média por corte) da aduela é, de acordo 
com a análise geométrica das deformações atrás 
realizada, 


Esta distorção está relacionada com forças de 
corte N., e N,, (Fig. 3), porque, como se verá 
adiante, estas forças são, cada uma, integrais de 
tensões tangenciais proporcionais à distorção 7,. 
Por sua vez o deslocamento U é igual ao 
deslocamento tangencial dos arcos devido as 
várias cargas actuando nelas mais a força tangen- 
cial horizontal 


N.. 
O Rai o (16) 


E o deslocamento V é igual ao deslocamento 
vertical da consola devido às várias cargas 
incluindo as forças verticais 


noam d Nxy 
p= (17) 


Poder-se-ia admitir, como foi feito noutros 
métodos ', que t, e t, eram independentes e 
então tentar obtê-los como forças separadas 
internas iguais, mas opostas, nos arcos e nas 
consolas, por meio de ajustamentos diferentes 
dos movimentos verticais e tangenciais. Este 
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procedimento, todavia, não é correcto porque 
ele despreza as relações fundamentais entre as 
forças de corte, as distorções por corte e as 
derivadas dos deslocamentos U e V (Eg. 15). 

Porém as duas forças de corte N,, e Nyx € OS 
dois momentos de torção M,, e M,, não são 
independentes. Como eles são os integrais das 
tensões tangenciais nas faces horizontais e verticais 
da aduela, basta apenas considerar uma das 
forças de corte e um dos momentos de torção 
para determinar os outros dois. Assim, admitindo 
que as tensões tangenciais seguem uma lei de 
distribuição linear, ao longo da espessura da 
barragem (Fig. 4), e isto sem perda de generali- 
dade uma vez que se poderia escolher outra 


Fig. 4 — Equilíbrio das tensões tangenciais 


distribuição qualquer, arbitrária, chamando 7,, 
T, € 7, respectivamente aos valores da tensão 
tangencial a montante, a jusante e no centro da 
secção, resulta que a tensão tangencial 7,, no 
ponto à distância z deste centro, é: 


T + To— T 
s u s u 
O 


x 2 e (18) 


Então nas faces verticais da aduela, resultará 


+ — 
E det (19 
Nç= | fada ma great, 0 
a: À 
2 
Ê (20) 
My= [ 2d (1,1) 


e 


2 
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e nas faces horizontais, sendo R o raio de curva- 
tura horizontal da superfície média da casca, será: 


Nu=) -E)de= 


o. (21) 
ts + Tu e 
4 “E 

M [ - (1 Z d E 
yZz z — R) Zz aZ — 

e (22) 

e, oh e 
e Hd To A s  53R 
A partir destas expressões, fazendo 
T=Ny 
e M, ai M ,x (23) 


2 2 
e eT 
Mo=My + New 7k Metqgr (24) 
e 
ê M 
N4=Ny—-—E=T-S a (25) 
R 12R R 


As eg. 14, 16 e 17 tomam então a seguinte 
forma : 


forças tangenciais (nos arcos) 


9N 
t=—E Ea (1 ps 
dy dy 12R“ R |, 
forças de flexão (de eixo vertical) 


M, 
x dy 


m (27) 
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ou se se introduzir a parcela devida a excentri- 
cidade de N,, (Eg. 10) 


9M, 
m',=-— + Ti (28) 
dy 
forças verticais (nas consolas) 
JT 
= ne 29 
E a (29) 
forças de flexão (de eixo horizontal) 
M , à Í eT 
y 
m=——— = M, + — (30) 
é dx dx 12R , 


No caso da barragem apresentar espessuras 
muito pequenas em comparação com o raio dos 


e 
arcos, os termos que contêm as relações R 
1 
e Rr podem ser desprezados e então 
xy = M x po M, (31) 
Na Ne = I 


e as Eg. 26, 27, 29 e 30 tomam a seguinte forma 
mais simplificada 


o T dT 
— e + p (32) 
dy d X 
dM, 9M, 
= € me (33) 
dy d x 


Note-se que esta simplificação não é indispen- 
sável para aplicação deste método e só tem sido 
usada naqueles casos em que a ordem de gran- 


[ss Rá é 
deza da relação E a torna legítima. Adiante 
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serão apresentados os erros cometidos quando 
se pratica esta aproximação. 


6 — CONDIÇÕES DE COMPATIBILIDADE 


O problema do cálculo das cascas, como 
se sabe “112 não se resolve apenas com a satis- 
fação das equações de equilíbrio e das condições 
limites ou de contorno, a menos que a casca 
seja encarada como membrana, desprezando-se 
os esforços transversos radiais e os momentos 
de flexão e torção, simplificação que já não 
é legítima nas barragens, por mais delgadadas 
que sejam. Torna-se por isso necessário juntar 
às equações de equilíbrio outras equações de com- 
patibilidade recorrendo à análise do comporta- 
mento elástico da estrutura. Se nas cinco primei- 
ras equações de equilíbrio atrás apresentadas 
(Eg. 9) se substituir N,, eM,, respectivamente 
por TeM,eN,, e M,, pelos segundos mem- 
bros das Eq. 24 e 25 passamos a dispor de cinco 
equações de equilíbrio envolvendo as oito incógni- 
tas seguintes: N,, Q,, No: Q,. M ,, M, jTeM,. 
Para a resolução do problema é necessário, por- 
tanto, introduzir mais três equações que relacio- 
nem deformações e tensões ou deslocamentos 
e forças. Essas três equações deverão impor 
a igualdade de todos os deslocamentos lineares 
e rotações sofridas pelas aduelas quer quando 
encaradas como pertencendo a arcos quer a con- 
solas. As três novas equações de compatibilidade 
a introduzir para levantar a indeterminação hiper- 
estática exprimem: (I) a igualdade dos desloca- 
mentos radiais; (Il) a igualdade das rotações 
de eixo vertical; (III) a igualdade expressa pela 
Eg. 15. Recordando a análise feita no cap. 3 
e admitindo como válida a aproximação resul- 
tante do facto de se desprezar as derivadas das 
distorções angulares 7, e d,, devidas aos cortes 
transversos radiais Q, e Q, concui-se que as 
igualdades (1) e (II) garantem a compatibilidade 
dos deslocamentos radiais W e das rotações... 
Yav' Yyx € “xy: A igualdade (III) garante a compa- 
tibilidade dos deslocamentos U e V e das rota- 
ções v,, e “,,- Este conjunto de oito equações 
satisfaz, portanto, integralmente, as condições de 
equilíbrio das forças e de compatibilidade dos 
deslocamentos. Note-se que a consideração das 
derivadas das distorções angulares é, e 2, (Eg.7) 
não introduziria nenhuma nova incógnita por 
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que são exclusiva consequência dos esforços 
0,eQ0.. 

Retomando agora as seis forças fictícias q,, Jc, 
Cita, € m, (Eg. 10 e 12) observa-se que 
elas, em conjunto com T e M, são suficientes 
para a determinação dos esforços N,, N, OQ 
O, M, e M,. Por outro lado as forças TeM, 
dão a conhecer t,, t y mem, por intermédio 
das Eg. 26, 28, 29 e 30 e q, é conhecido a par- 
tir de q, (Eg. 12). Então o problema reduz-se à 
determinação de gq, TeM,. 

Estas seis forças fictícias permitem a determi- 
nação de deslocamentos da barragem. As três 
das Eq. 10, nos arcos, produzem deslocamentos 
radiais W,, deslocamentos tangenciais U, e rota- 
ções O, de eixo Vertical, de acordo com as fór- 
mulas bem conhecidas da teoria dos arcos. Não 
é necessário calcular os deslocamentos verticais 
e as rotações de eixo horizontal dos arcos. Os 
deslocamentos verticais dos arcos são tomados em 


7) 
consideração no cálculo através de — = Vos 
x 


que é considerado como uma variável indepen- 
dente, e uma vez que se imponha que as rotações 
de eixo vertical dos arcos sejam iguais às rotações 
de eixo vertical das consolas, a compatibilidade 
das rotações de eixo horizontal torna-se redun- 
dante. 

As outras três forças fictícias correspondentes 
às consolas permitem a determinação dos deslo- 
camentos radiais W., dos deslocamentos verti- 
cais V, todos funções de q, t, em, e as rotações 
de eixo vertical 9, (devidas apenas ao M,). Os 
deslocamentos tangenciais horizontais das conso- 
las não são calculados mas não tomados em con- 
po = OU 

ue dy Í 
considerada como uma variável independente. 
Como já se referiu, as as forças fictícias m, Mm, 
t, e t, são derivadas de funções em que intervêm 
as incógnitas M, e Tea somade v,, e v, é 
proporcional a T, pelo que há apenas três variá- 
veis a determinar que são g,, M, e T. Todavia, 
por razões de simplicidade e rigor , tomam-se, 
normalmente, as seguintes quatro variáveis: 9,, 
M, Te 4- 

Também por razões de simplicidade e rigor 


sideração através da derivada 


a condição de compatibilidade (111), y, = = 
e 
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ot 
= 2U + AN exprime-se pelas Eq. 34 e 35 que 
dy dx 


se discutirão a seguir. 


A compatibilidade dos deslocamentos verticais 
é satisfeita pela expressão: 


OdV E dV 
v=— oq =p =D (34 
ei dx Ge Yyz dx 8a) 


sendo V o deslocamento vertical das consolas. 
Quando se analisam os deslocamentos tangen- 
ciais horizontais, pode-se notar que os centros de 
gravidade das aduelas, que constituem uma con- 
sola, formam uma curva torsa no espaço, 
A direcção do eixo tangencial x varia, por con- 
A * . 
sequência, ao longo da linha de centros e a 
dU 
m- 4 RS = A . o 
expressão v,, = —— é apenas válida na vizi- 
0y 
nhança do centro da aduela. Se a curvatura ver- 
tical e as variações de direcção do eixo xx são 
consideradas, a equação de compatibilidade dos 
deslocamentos tangenciais, será expressa pela 
equação: 


dU Ui wºt 


— 35 
dy R dy ama 


J geme 
T yz 


sendo R o raio da superfície média e 6 o ângulo 
horizontal entre o eixo da aduela e o eixo da 
barragem. 

Deste modo, podem escrever-se as seguintes 
equações de compatibilidade : 


Ve, (Q ui mo.) e W. (Q-Ja P, by, M,, m) 


0 (da te mo) = 0. (M) 


a + 
U, (Jar Le m',) + fu, (Ja le m,) E = 


) A 
us ds W, (Ja Er m,) e Ds | dy=0 
dy 


T 
A (q — Ja P, ty m, m,) — [tt dx == O 


A ça TED, (36) 


dy 


x 


341 


d eT 
m, — — + ——) = 0 
' q IM, 12R' 
z 

e! a a 
dy 12R R 
91 

a NR 

y d x 


A solução deste sistema permite a determina- 
ção de oito incógnitas em cada aduela q, M,, 
TV» MD, t Et É este sistema de equa- 
ções que resolve o problema da determinação 
das forças internas das aduelas. Normalmente, 
para o cálculo de uma barragem simétrica, 
usando uma malha de seis arcos e seis consolas, 
o número de equações é de 189. 

Calcula-se em primeiro lugar o efeito das 
forças Ga, GC (= q— qu), My ty, tys Me m, nos 
arcos ou nas consolas para cargas unitárias, e 
depois introduzem-se esses efeitos em termos de 
coeficientes de influência k nas aduelas de cruza- 
mento de arcos e de consolas, como se indica 
em seguida. Nas expressões que relacionam 
tu ty mo em, comM, e T, substituem-se as 
derivadas por diferenças finitas. 


(Wa) = : (kJ (qa) + (8) (mo) + 
+ : (k), (t,); 


(U,), = idêntica 
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(37) 


(04), = idêntica 
(No), = 2 (k), (g— qu); + a (k), (m,), + 


-t- 2 (k),; (t,); 
(V.), = idêntica 


(96) mm : (e), (Mp; 


Deste modo, as consolas necessitam de ser 
analisadas para séries de forças triangulares 
unitárias: radiais, q.;, verticais tyié flexões de 
eixo horizontal, m,,, e torções verticais M,,. Não 
é necessário considerar forças tangenciais. 
Semelhantemente, os arcos devem ser analisados 
para cargas radiais q,;, tangenciais t,, e flexões 
de eixo vertical m,,, e as cargas de torção de 
eixo horizontal não necessitam ser consideradas 
nos arcos. 

A partir da solução do sistema de Eq. 36 
determinam-se os valores de cada uma das cargas 
Gs ta ba MD. m,, Me T em cada aduela de 
cruzamento. A partir destas cargas, determinam-se 
as forças internas que entram nas Eg. 11 e 13. 
Pode então fazer-se a determinação completa das 
tensões nas faces e no interior da barragem e 
também das reacções da barragem contra a 
fundação. As fórmulas usadas para a determi- 
nação das tensões principais são as que se 
deduzem fácilmente das condições de equilíbrio 
de tetraedros junto aos paramentos. 


(Continua) 
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LINHAS POLIFÁSICAS UNIFORMES E LINEARES EM FUN- 
CIONAMENTO SINUSOIDAL PERMANENTE EQUILIBRADO 


II— SIMETRIA, PASSIVIDADE E RECIPROCIDADE 


RESUMO 


Apresentam-se as expressões gerais de uma linha poli- 
fásica uniforme finita e definem-se as condições de sime- 
tria, passividade e reciprocidade, discutindo-se as res- 
pectivas consequências nas grandezas fundamentais 
da linha, 


1 — INTRODUÇÃO 


O estudo geral das linhas uniformes de trans- 
missão por meio do cálculo matricial conduz 
a expressões que convém analisar no caso prático 
das linhas finitas e quando se verificam condições 
especiais de passividade e reciprocidade, pois 
nestas circunstâncias as propriedades gerais 
assumem aspectos simplificados. 

Na presente análise considera-se fundamental 
o conhecimento das leis gerais aplicáveis às 
linhas polifásicas uniformes e lineares (1). 


2 — LINHA UNIFORME FINITA 


2.1 — Equações da linha finita 


Seja 1 o comprimento de um troço de linha em 
que os valores das tensões e das correntes no 
terminal de entrada (extremo ligado ao gerador 
polifásico que a alimenta), isto é para x==o, 
são V,=V, e l,=1,e os valores correspondentes 
no terminal de saída (extremo ligado à carga ou 
receptor), quer dizer, para x =1, são V=v, e 
|=1,. 

Pelas equações hiperbólicas matriciais da linha 
podem exprimir-se as grandezas à entrada da 
linha em função das homólogas à saída (sistema 
(43) em [1]) 
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The essential expressions of a finite homogeneous 
bolyphase electrical line are presented, The phisical and 
mathematical conditions of symmetrical, passive and 
reciprocal lines are defined and the main caracteristics 
are discussed. 


V ch G (sh G)Z, o” 
(|| 2;'n6 2z-'(h6)z]| 1, | 
Define-se então a matriz de tran. jerência da linha 


finita 
— [eh(T) sh(TDZ, 
FE pimiry Z-atriz| 2 


visto que a matriz de propagação das tensões tem 
o valor 


GT. (3) 
A equação geral da linha finita será 


N=TN, (4) 


onde N, representa o vector tensões-correntes 
à entrada da linha, N, é o vectór tensões-corren- 
tes à saída da linha e T traduz a matriz de 
transferência do troço de linha de comprimento 
|, ao qual correspondem as condições de saída N, 
quando na entrada as condições são N,. 

Vê-se assim o significado da matriz dé transferência 
da linha finita: matriz que relaciona as grande- 
zas de entrada com as grandezas de saída de um 
troço de linha uniforme, por intermédio das 
matrizes das impedâncias de onda das tensões e 
das constantes de transmissão das tensões. 
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2.2 — Quadri-multipolo equivalente da linha 
finita 


Considerando a linha finita como um quadri- 
-multipolo, constituído por 4m terminais, sendo 
2m de entrada e 2m de saída, equivalente à 
linha polifásica, a equação da linha pode tomar 
o aspecto híbrido 


relativo às equações 


V,=AV,+8BI, (6) 
| =CV,+0DI, (7) 


que definem a matriz de transferência 


“[aAgSB 
=| | (8) 


com os valores determinados por 


A=ch G (9) 
= (sh6)Z, (10) 
=Z, (sh6) (11) 
D=2"' (ch 6)Z (12) 
ou 
p= ABI ME Es 05) 
2! 4! 


B =|! à a PZva E YA... fizas 


c=[1 4 e PYZ+ = (YZ) + e [YA 


! 


Ds PV A+ (UI 4. (16) 


Por razões equivalentes a matriz de transfe- 
rência inversa será 


“A -B 
h=| E + (17) 
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Análogamente se pode raciocinar com a matriz 
de transferência T' referida às grandezas das 
correntes definindo-se então as constantes gene- 
ralizadas referidas às correntes À”, B', C' e D'. 


2.3 — Impedâncias de onda 


De um modo geral, para se determinar a impe- 
dância de onda de uma linha finita não longitudi- 
nalmente simétrica, fecha-se a linha por uma 
rede de cargas com 2 m terminais de matriz 
de impedâncias Z, e mede-se à entrada a matriz 
de impedâncias Z,; depois fecha-se a linha 
à entrada por uma rede de cargas de matriz 
de impedâncias Z, e mede-se na extremidade 
de saída uma matriz de impedâncias Z,. Estas 
duas matrizes de impedâncias, que satisfazem 
a estas propriedades, são as matrizes de impe- 
dâncias de onda respectivamente à entrada e à saída 
da linha finita. 

Num sistema físico com simetria longitudinal 
(como é o caso da linha uniforme, como se verá) 
pode simplificar-se a determinação da impedância 
de onda: supondo a linha finita fechada por uma 
carga de 2m terminais, ligados cada um a um 
terminal de saída da linha, a matriz das impedân.. 
cias de onda é caracterizada pelas impedâncias 
de carga tais que a matriz das impedâncias da 
linha e da carga aos seus terminais é igual 
à matriz das impedâncias de carga. 

Deste modo conclui-se que V,=Z,l, e 
V,=Z,],. Nestas condições quando a linha 
finita está fechada pela sua impedância de onda 
anulam-se as parcelas relativas às ondas de 
reflexão e portanto 


V,=eºyv, (18) 


e a Me 2 MM (19) 
Outro modo de definir as impedâncias de onda 
consiste em calcular a matriz das impedâncias 
à entrada da linha infinita, como se verá. 
No funcionamento da linha finita em vazio, 
caracterizado pela condição |, == 0, as equações 
contidas em (1) reduzem-se a 


V,=ch6.v, 
e =2Z, 'hG.v, 


(20) 


(21) 
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que condensam, respectivamente, as interpreta- 
ções físicas das constantes generalizadas À e C. 


De (20) e (21) 


V,=(cth 6) Z,.l, (22) 
donde a matriz das impedâncias à entrada da 
linha em vazio 


Z,, = (cth 6) Z,. (23) 

No funcionamento da linha finita em curto-circuito, 

condicionado por V, = 0, analogamente se pode 
verificar 

Z, = (th 6) Z,. (24) 

Utilizando as equações referidas às intensida- 
des obtém-se similarmente 


Z, = Z, (cth 6º) (25) 


2 =2, 6h 64, (26) 
A analogia destas expressões sugere que haja 


uma relação entre Z ze Z : . De facto de (23) e 
(25) 


(cth 6) Z, = Z, (cth 6º) (27) 
e de (24) e (26) 
(th 6)Z,=2Z, (th 6º) (28) 
pelo que 
Z,=2Z, (29) 


e, substituindo em qualquer das expressões (27) 
ou (28), 


(30) 
que, pelas respectivas definições, conduz a 
p pras; [ 


(31) 


Deste modo a matriz das impedâncias de onda 
pode-se representar por qualquer das expressões 


| 
Zo=(ZY) *Z 2) 


Z.=(2Y) “yr! (33) 
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— E! 
ou, sucintamente, 
É ma. 
Z =7 2y 2 (36) 
o) 


2.4 — Simetria longitudinal 


A condição de simetria longitudinal da linha 
caracteriza-se pelo facto da permutação dos 
índices da equação geral da linha exigir única- 
mente a inversão da matriz de transferência, 
isto é, N=N,eN,=T,N,com T,T”! pelo que 
TI,=1, donde, substituindo pelas expressões 
(8) e (17), se obtém as condições 


A?-BC=1 (37) 
D—-CB=1 (38) 
AB=BDe,sedetBo,D=B-'AB (39) 
CA=DBe,sedetC+o,D=CAC”! (40) 


As duas últimas condições mostram que as 
matrizes À e D são semelhantes, visto que se pode 
obter uma da outra por uma transformação de 
semelhança. 

É fácil mostrar que, dada a existência de sime- 
tria longitudinal, o determinante da matriz 
de transferência tem o valor +- 1. De facto 


T=| À ADE | 


Ê AR. Mt) 


pode decompor-se em 


SI E a 

b (A—1) C | E! (A—1) C [uz 
0 ] C IJ ] 

em que cada matriz tem o determinante igual 

a + 1, sendo portanto 


det T =+ 1. (43) 
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2.5 — Linha com simetria transversal 


A condição matemática de simetria circular 
exprime que as matrizes Z e Y são transforma- 
das em matrizes diagonais pela mesma transfor- 
mação diagonal. 

Como consequência facilmente se deduz que 


ZY =Y2Z 
P = 
G — 6" 
=, 
À = 


3 — LINHA PASSIVA 
3.1 — Condições de passividade 


A condição física de passividade de uma linha pode 
enunciar-se dizendo que a potência activa à entrada 
de uma linha passiva não pode ser negativa; 
pode ser nula se, estando a linha ligada às fontes, 
não houver dissipação de energia, isto é, se 
com 150 for P=o. 

Ora a potência activa é a parte real da potên- 
cia aparente complexa, P==Re 5, isto é, 


onde Vel, | São os valores da tensão e da 
corrente complexas à entrada da fase j, ou, em 
notação matricial, 


P=Rev l (44) 


ou seja P = Re Pp V, (45) 


e, sendo Z, a matriz das impedâncias à entrada 


da linha, 


P=Rel, Z,|, (46) 


quer dizer, P é a forma hermítica de S =1,, Z,|, 
e atendendo a que Z,=Z,, +Z,, isto é, 


Zu = (2 + Z, )/2, vem 
P=h Zmh- (47) 
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Analogamente, como a potência aparente 
S=V,,Y,,V, deduz-se uma expressão equiva- 
lente para a potência activa em função da matriz 
das admitâncias à entrada da linha, sendo por- 
tanto análogas as condições matemáticas para as 
formas hermíticas Z,y e Y,, para que corres- 
pondam a uma rede passiva. 

A condição física de passividade exige que 
P>0 ou P>0. Assim, exceptuando o caso em 
que |,=0 para o qual toda a forma hermítica é 
idênticamente nula, será (na hipótese |, + 0) 


IH 


P>o : matriz definida positiva 
P > O : matriz semi-definida positiva 
PO : matriz definida negativa 
P<.O |: matriz semi-definida negativa 
—c<P<<-+oo : matriz indefinida. 


Então a condição matemática, necessária mas não 
suficiente, de passividade exige que a parte hermítica 
da matriz das impedâncias (ou das admitâncias) 
à entrada da linha seja definida positiva ou semi- 
-definida positiva. 

Por outro lado, é condição necessária e sufi- 
ciente para que uma matriz hermítica seja definida 
positiva que todos os seus valores próprios sejam 
reais positivos e para que seja semi-definida posi- 
tiva que alguns valores próprios sejam positivos 
e outros nulos. 

O critério estabelecido é no entranto mais 
fácil de verificar através de certos determinantes: 
a condição necessária e suficiente para que uma 
matriz hermítica seja definida positiva é que todos 
os determinantes principais que se podem extrair 
da matriz sejam positivos (teorema de Syl- 
vester). 

Mas uma matriz hermítica semi-definida posi- 
tiva é singular e portanto os seus determinantes 
principais não são todos positivos, definindo-se 
então o critério seguinte: a condição necessária 
e suficiente para que uma matriz seja semi-defi- 
nida positiva é que, após permutação eventual 
das linhas e das colunas, os determinantes prin- 
cipais de ordem crescente sejam positivos até à 
ordem r, sendo todos os outros nulos, dizendo-se 
nestas circunstâncias que a matriz é semi-definida 
positiva de ordem r. 
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3.2 — Matriz das constantes de propaçação 


Na linha passiva desaparece a ambiguidade do 
sinal das matrizes ! e |'* introduzida pela extrac- 
ção da raiz quadrada na sua definição (1). Realmente 
os sinais devem ser tais que Z, Seja a matriz 
de uma linha passiva. Para isto é necessário que 
os valores próprios de [ correspondam aos únicos 
amortecimentos e desfasagens fisicamente pos- 
síveis. 

Dados Z e Y calcula-se a matriz ["=ZY e 
os valores próprios da matriz I'* dados por 

det (Tº—))=0 (48) 
que define um polinómio de grau m em À e por- 
tanto a equação característica cujas m raízes são 
os valores próprios da matriz I'”. Utilizando uma 
transformação de semelhança da matriz ['* por 
meio da matriz L obtém-se a matriz À semelhante 


: o. 
à matriz E 


LT! TP? L=A (49) 


onde À é uma matriz diagonal. Então = LAL! 
e portanto 


(50) 


sendo A'* uma matriz diagonal cujos elementos 
da diagonal principal são as raízes quadradas 
dos valores próprios de [”, isto é, + VA com 
k == 1,2,...,m,e os restantes nulos. A escolha 
dos sinais deve corresponder ao amortecimento 
crescente com x para todas as funções próprias, 
sendo portanto os valores positivos. 


3.3 — Matriz das impedâncias de onda 


Como se referiu pode definir-se a matriz das 
impedâncias de onda como a matriz das impe- 
dâncias à entrada de uma linha infinita, 

Sejam n linhas finitas em cadeia, com a matriz 
das impedâncias de onda das tensões Z e a matriz 
das constantes de transmissão G. De (23) deduz-se 
que a matriz das impedâncias à entrada da cas- 
cata de linhas iguais em vazio será Z,, = (cthnG) 
Z ,, e em curto-circuito Z,, = (thn G) Z,. 
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Fazendo tender o número n de células para 
infinito a matriz de impedâncias à entrada de 
uma linha uniforme infinita é, em vazio, 


Z voo = lim (cth G) Z, (51) 
1 —> CO 
e em curto-circuito 
Z, «co =lim (thn 6) Z, . (52) 


—» 09 


No caso de existência de passividade as matri- 
zes Z e Y satisfazem às condições de matriz 
definida positiva, assim como a matriz Z,. Nestas 
circunstâncias verifica-se a convergência das rela- 
ções anteriores pois 


lim (cthn G) =lim (th n G) = 1 


—» CO 


(53) 


 —> 00 


como se pode verificar, por exemplo, pela relação 
thn G=(I—(chn G)-2) !!2, Em conclusão 
DG dad E pi aa (54) 

que exprime o significado da matriz das impe- 
dâncias de onda como sendo a matriz das 
impedâncias à entrada de uma linha infinita 
passiva. 

Note-se que pode haver casos em que a matriz 
das impedâncias de onda indica o caracter pas- 
sivo da linha sem que existam os limites (53). 


4 — LINHA COM RECIPROCIDADE 
4.1 — Condição de Reciprocidade 


No caso de uma linha exclusivamente consti- 
tuída por dipolos (resistência, capacidade e auto- 
-indução) e por transformadores com relações de 
transformação reais, as matrizes impedância e 
admitância aos terminais da linha são simétricas. 

A reciprocidade de uma linha traduz-se mate- 
maticamente pela simetria das matrizes impe- 
dância e admitância aos terminais, Z = Z, e 
Y = Y,, tendo como consequência a simetria das 
matrizes das impedâncias unitárias e das admitân- 
cias unitárias dada a simetria longitudinal. 

O princípio de reciprocidade só se aplica a 
linhas em que as perdas joule, a energia magné- 
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tica e a energia electrostática sejam funções de 
estado para todos os valores instantâneos das 
tensões e das correntes aos terminais | 2 |. 


4.2 — Impedâncias de onda 


Como Z e Y são matrizes simétricas verifica-se 


ZY=(Y 2), (55) 
e portanto 
Ls (56) 


quer dizer, as matrizes das impedâncias de onda 
das tensões e das correntes não simétricas. 

Por outro lado, a simetria de Z e de Y implica 
a simetria de Z,, pois 


Z=2 (YZ) "= 
Es FAR o DO it 


2 t? y)! em 
=((2Y) “Z), = 
= Zw- (57) 


Este resultado era de prever conjugando (56) 
com Z,= Z,. 

Note-se que a simetria de Z e Y não implica 
a comutatividade Z Y ==Y Z pois o produto de 
duas matrizes simétricas não é em geral uma 
matriz simétrica. Mas os produtos de duas 
matrizes simétricas têm a mesma equação carac- 
terística e portanto os mesmos valores próprios, 
propriedade muito importante para simplificação 
dos cálculos numéricos. 


4.3 — Constantes de propagação 


Quando as matrizes eZY são simétricas 
a matriz das constantes de propagação das cor- 
rentes é a transposta da matriz das constantes 
de propagação das tensões e portanto também 


são simétricas 
r=(y 7)º= 
am (Y, Z)!'2 =— 


= (20), = 
= K, (58) 


e por raciocínio análogo 


G = 6, (59) 


348 


4.4 — Passividade 


Numa linha com reciprocidade a parte hermé- 
tica da impedância aos terminais coincide com 
a sua parte real pois a transposição não as modifica 


Zw=> (2 +21) =ReZ,=R (60) 


e portanto a passividade da linha será determi- 
nada pelos valores próprios de R. 


5 — CONCLUSÕES 


A teoria das linhas polifásicas revela que uma 
linha uniforme e linear apresenta simetria longitu- 
dinal e goza das seguintes propriedades: 


a) As matrizes das impedâncias de onda das 
tensões e das correntes são iguais (Z, — Z,); 


b) As matrizes das constantes de transmis- 
são das tensões e das intensidades são iguais 


(G = G"); 


c) As matrizes das constantes de propaga- 
ção das tensões e das intensidades são iguais 
(=; 


d) As matrizes de transferência referida às 
grandezas das tensões e das correntes são iguais 


rss Te): 


A passividade conduz à univocidade de cálculo 
da matriz das constantes de propagação e per- 
mite definir a matriz das impedâncias de onda 
como a matriz das impedâncias à entrada de uma 
linha infinita. 


A reciprocidade implica propriedades de simpli- 
ficação no estudo e cálculo quer das impedâncias 
de onda quer das constantes de propagação. 
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COMPUTER CALCULATION OF STRESSES 
AND ROCK FAILURE ANALYSIS 
IN SOME BLASTING SITUATIONS 


RE SUMO 


Vários programas de computação são apresentados 
para o cálculo de tensões indusidas no interior de rochas 
submetidas à detonação de cargas explosivas esféricas 
e cilindricas. Analisam-se diversas geometrias e as dis- 
tribuições de tensões são obtidas em função do tempo. 
À energia de deformação associada à onda de choque 
emitida é considerada e apresenta-se também um pro- 
grama de computação para a mesma, A análise da rotura 
da rocha é estabelecida com base num critério dinâmico 
digno de confiança. Os resultados obtidos estão em 
concordância com dados conhecidos de arranque de rochas 
em trabalhos mineiros e experimentais, 


1 — INTRODUCTION 


The stresses induced in the rock surrounding 
a concentrated charge of explosive, after its 
detonation, can be determined by means of the 
mathematical theory of elasticity with good 
approximation, because in such cases the devia- 
tions from perfectly elastic behavior are small. 
The usual approximations of elastic wave theory 
are valid for most materials and most problems. 

This type of stress analysis was first deduced 
for the case of the explosion of a confined 
charge within a rock mass (!) and latter for the 
same phenomenon, but in the vicinity of one 
free face (*). 

Therefore, a computer program was reported (*) 
in view of the complexity and mass of calculations 
necessary to accomplish the stress determination 
for the second of those cases. 

More recently, a very important hypothesis 
have been introduced in the theory of blasting (*) 
by means of which a cylindrical charge of explo- 
sive is assumed to be equivalent in effect to a 
composite of concentrated charges. Essentially, 
the stress pulse from a cylindrical charge is 
simulated by the superposition of elemental 
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SYNOPSIS 


Various computer programs are reported for the 
calculation of stresses within a rock mass subjected to 
the detonation of spherical and cylindrical charges of 
explosive. Different geometries are analysed and stress 
distributions are obtained as transient solutions. Refe- 
rence is made to the radial strain energv associated to 
the stress wave coming out of the detonation and 
a computer program is also presented for its calculation. 

Rock failure analysis is deduced in terms of a relia- 
ble dynamic failure criterion. Results agree with known 
data from mining and experimental blasting operations. 


stress pulses coming from a line of concentrated 
charges that make up the cylindrical charge. 
This assumption, in which theoretical treatment 
has produced good results in agreement with 
experimental field work (”) provides a good basis 
for the mathematical treatment of the problem, 
suggesting its application in new methods of 
design blasting patterns, featuring column charges, 
by far the most employed procedure in practice. 

This paper presents the extension of the 
referred computer program to the case of a 
cylindrical charge detonating in the vicinity of 
two free faces of rock, (the so called bench blas- 
ting), assuming that composition of effects. The 
transition is achieved through an intermediate 
step, describing the case of a spherical charge 
detonating near two free faces. This phase was 
developed more in connection to the assessment 
of the mathematical model of wave reflection, 
rather than some possible practical applications. 

In all the cases, the purpose of the method 
consists in establishing the induced stress distri- 
bution (represented by the magnitudes and 
directions of the principal stresses) over small 
increments of time. Then, the correlation with 
an appropriate rock failure criterion (maximum 
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dynamic tensile strength of the material) is set 
up for the prediction of the produced fractures 
and crater boundaries. 

Results show good agreement with known 
values from practical cases of mining and both 
the laboratory and field experiments reported in 
the bibliography. 

In particular, we emphasize the enormous 
capability of the computer in order to obtain 
the different stress distributions when slight 
changes in the input data are introduced. This 
procedure is extremely convenient to investigate 
the separate influence of the various parameters 
in the phenomenon, indicating clear trends for 
studies of optimization, as well as supplying 
valuable information for the explanation of 
mechanisms of rock fracture, as transient situa. 
tions. 

Therefore, the method can provide an impor- 
tant tool for the scientific design of blasting 
operations. 


2 — DETONATION OF A SPHERICAL CHARGE 
IN AN INFINITE MASS OF ROCK 


Considering a homogenecus, elastic and infinite 
mass of rock within which detonates a spherical 
charge of explosive, the principal resultant stress 
wave is a longitudinal wave propagated through 
the surrounding medium with a velocity C, cha- 
racteristic of the rock. 

Using polar coordinates centered at the charge 
center E, each point A in the rock mass is subjected 
to a certain displacement U, with the components 


(Fig. 1): 


T represents the time from the instant of deto- 
nation and r the distance EA. 

The theory of elasticity shows that the stresses 
at the point A are defined by 


3U e é 
e. =(A +24) — + 24 — 
JT r 

R$) . U 

G() o Che À — + 2(4 + 4) Re 

dr r 


Cn=€, = t,.=0 


tr 


350 


where ) and vu are the Lamé's constants of the 
rock. 


| 
| 
I 
| 
sd 


Fig. 1 — Polar coordinates of a point A in the interior 
ofthe rock mass. 


The displacement changes with time can be 
expressed as the product 


U,=U (r,T) =U, (r). Uw (T) 


where U, (r) represents the peak value of the 
particle displacement, which is a function of the 
distance r, and Uw(T) is the change in displa- 
cement with time, commonly known as the wave 
shape for displacement. 

The arrival of the longitudinal wave to the 
point A will take place in the instante T,, 
(i stands for incident and p for P-wave) where 


and the corresponding stresses are, using the 
same notation 


9Ul(r,T ) Vig E 
Cio fd dp) ——— P ) ip) 
dT r 
UC à Ego) HtE 
op up á + 2(4+ uv) gi 
Crcip = Cotip = Cp =0 


The derivative terms can be written: 


dUl(r,T) ã 
aU, ( 9 Uw(T, 
ma À wir) E qm 
r 
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and the last term can be transformed because 
O Uw( Ti) dUw( Ti) O O To 
dr EF T, or 


dUwíTo) dT, 


/ = à 
(Tc) 


dT, OT 
0 1 dUwíT;) 
Cp d Ts 


Introducing the radial particle velocity 


OU (r, Tp) 
oT 


V(r, Ti) = 


and substituting, we obtain 


9U(r,T dU Vir, 1...) 
GT) 2 O ar) To) 


dr dr ô 


Thus, the stresses at the point A, such that 
r == EA, are: 


p= 0+2u | É Uva) — 
á vem 123 VET) 
e, 1 r 
dU. (r) 
*qip = Sp = À | ri Uw (Ip) — 
V tr, 1 U (r,1,,) 
A +20 + 0) 
e, | r 
Croip — Cotip = Crrip o 


Because of the absence of shear stresses, 7 
7 sip and 2, are principal stresses. 


Hp” 
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I is then necessary to define the quan- 
tities U,, Uw and V, including their variation 
with time and distance. 

The attenuations of maximum displacement 
and maximum particle velocity are known 
to follow inverse power functions of the distance. 
Theoretically, in a perfect homogeneous medium 
the decay of stress waves is found to be exponen- 
tial, with a coefficient dependent on the internal 
friction of the solid. 

In blasting, these exponents (m and n) have 
larger values in the vicinity of the shot point, 
but with reasonable approximation the following 
laws are valid. (*) 


Utr, Tp) = A, 1 " Uw (To) 


V(r, Tp) =B, ro "Vw(T;) 


where A, and B, are constants (intercepts) and 
Vw (T,,) is the change in particle velocity with 
time (wave shape). 

It is evident that 


— 
U, (1) = Agr 
and 


V()=B,r;” 


and these are the necessary relationships for the 
stress calculations, providing that the experi- 
mental measurements give reliable values for the 
constants involved. 

A good description of the techniques and 
results obtained for this purpose is found in (3) 
and it makes possible to accomplish the nume- 
rical calculations. A summary of the field 
conditions and characteristic parameters is pre- 
sented next. 

The rock is a ferrous ore (mainly composed 
by magnetite and specularite) with an average 
specific gravity of 3.8 g/cm* and the explosive 
selected is Geogel 60º/, characterized by a deto- 
nation velocity of 6.3 >< 10º cm'sec. 

Time dependent displacement and particle 
velocity in that roch were obtained from field 
array experiments. Table No 1 gives the prin- 
cipal data needed for the calculations. 
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TABLE No 1 


Variable Fortran Value Description 
Symbol 
Cp CP 64> 10º cm's Velocity of propagation of the longitudinal waves 
Cs cs 3.42>X 10º cms Velocity of propagation of the shear waves 
/ AL 6.25>10H dyne/cm?| First Lamé parameter for the rock 
a AM 4.15>< 10!! dyne/cm? | Second Lamé parameter or shear modulus of the rock 
A q AO 40.3 cm Displacement intercept (for Geogel 60"/,) 
Bo | BO 3.36 >< 10º cm/s Particle velocity intercept (for Georgel 60/,) 
n EN 1.4 Decay exponent of radial displacement in the longitudinal wave 
m EM 1.4 Decay exponent of particle velocity in the longitudinal wave 
h | EH A Decay exponent for displacement in the shear wave 
k EK 1. Decay exponent for particle velocity in the shear wave 


The wave shapes for displacement Uw and for 
particle velocity Vw characteristic of the stress 
wave induced by a shot of Geogel 60 º/o within 
that rock, are depicted in Fig. 2. These data were 
obtained (º) after the recording of oscillogram 
traces and revealed a small variation in the wave 
shapes with distance. Therefore, we assume that 
they represent with sufficient accuracy the 
corresponding functions over a wide range of 
distances from the shot point. 

According to this set of elements, a computer 
program was established (Program NORE) and 
it is presented in the Appendix 1. 


Time [IU"* sec) 


Fig. 2—- Wave shapes for displacement (Uw) and for 
particle velocity (Vw) used in the computer program. 


Various distances to the center of the detona- 
ting charge were considered in the rock for the 
sake of calculation of the principal stresses, over 
30 equal time increments, each one of 107* 
seconds. Tension is regarded as positive and 
compression as negative stress. 
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Fig. 3 contains the conclusions of the output 
data for induced stresses in function of time and 
at increasing distances from the shot point. 

The radius of the spherical cavity containing 
the explosive is 8.2 cm (3 1/4”) for this and the 
next blasting situations. 


def 


Time (msec! 


Q r= AE =50em 
e lc 
Tlel th 
2 x10º E. 
Hime imsec] 
Ú r = AE =WUcm 
r=ÃE-Z00cm 


é x10" 


Hime lmse 
f AE = 30bem 


Fig. 3— Stress distribution variation with time at diffe- 
rent distances from the explosive charge. 
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In each of the diagrams, the maximum dyna- 
mic tensile strength of the rock (200 x 10º 
dyne'cm?) is referred, according to (3). 

The extension of fractures with time is assumed 
to follow the path of sucessive points where the 
value of that strength is exceeded. The cracks 
follow the direction of o, when 9, is greater 
than the dynamic tensile strength, and they ter- 
minate when these values are equal. By means 
of this assumption, that provides good correla- 
tion with the results of blasting tests (7), (*)itis 
possible to predict the extension of the fractured 
zone of rock surrounding the shot, as well as 
the crater shape and volume, when free faces 
are in the proximity of it. 

An interesting correlation between the exten- 
sion of radial cracks and the strain energy asso- 
ciated with the stress wave coming out of the 
detonation, can also be deduced from the com- 
puter program. 

The radial strain energy W, radiated outward 
per unit area at the instant t is (º), ('): 


3 t 
= PC 
W,= É a! 2 da 
à o 


where q is the material density, C, the longi- 
tudinal propagation velocity and «, the radial 
strain induced, in this case given by: 


0U, 
E — <a 
: dr 
dU,(r V(r,T..) 
dU, (1) U Ti) = ip 
dr Cp 


For a certain distance r from the shot point, 
the total radial strain energy is: 


and if r represents the radíus of the spherical 
volume cracked around the shot (sometimes 
called equivalent cavity) we obtain the amount 
of strain energy required for the radial cracking 
mechanism produced by the stress wave. 

From the results of the program NORE and 
after the application of the tensile strength cri- 
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terium this radius is known to be approximately 
r==190cm, and at this point the corresponding 
instant when the maximum tensile stress indu- 
ced exceeds the dynamic tensile strength of the 
rock is t=7. 0>< 10-* seconds. 


These quantities and the values of p and C, 
mentioned obove yield for W, the amount of 
1.432>< 10!! erg'cm? and a total strain energy 
required to fracture the rock within the equiva- 
lent cavity, of W, = 6.491 10*º erg. 

The complete calculation of strain energy was 
also obtained for a wide range of distances from 
the charge and over 30 increments of time for 
each of the distances. The program is also 
presented in the Appendix 1 — Part II (program 
ENERGY). 

I can be assumed from the analysis of the 
output data that the constant decrease of strain 
energy with distance for the same instants, 
coupled with the empirical evidence of crack 
development, suggests the existence of a critical 
strain energy below which the fracturing stops. 

The variations strain energy-time for various 
distances from the shot point are presented in 
Fig. 4. The assumed maximum distance of frac- 
ture is plotted, and it separates the two regions 
of the rock mass. 


Wo | queseeeeseneneenemerenecencenterceranerrereransersacas r= 40 em 
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Fig. 4 — Strain energy variation with time and distance- 


The subcritical values of the strain energy are 
related with zones of the rock mass where the 
rock remains intact. Therefore, the two criteria 
of failure can provide, independently, the infor- 
mation needed for fracture prediction when their 
values are compared with some critical charac- 
teristic parameters of the rock. 

In particular, if the critical strain energy is 
known a priori, it could be used in a much more 
reliable way to predict crack extension, than the 
dynamic strength alone. 
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3 — DETONATION OF A SPHERICAL CHARGE 
IN THE VICINITY OF ONE FREE FACE 
(CRATERING) 


The method of stress analysis reported for 
the calculation of magnitude and direction of the 
principal stresses in the case of blasting with 
a spherical charge near a free face of rock, uses 
the plane wave reflection theory. Again, the 
deviations from the real case of spherical waves 
are not important in view of the relatively redu- 
ced distances of propagation involved in this 
case. 

The stress at any point A within the rock 
can be synthesized by transferring the stresses 
due to the incident longitudinal wave, the stress 
due to the reflected longitudinal wave and the 
ones due to the reflected shear wave. As these 
three waves arrived at A in distinct instants (and 
following that order), the respective composition 
is done with consideration of their correspondent 
times of arrival to the point. 

The situation is presented in Fig. 5. 


pes EA 
r,=ÃAÃO, 
,=AÃO, 
r,=0,x; 
=x, 


Elxe,yE) ly 


Fig. 5 — Polar coordinates system used for the case 
of blasting with one free face (free face is the 
plane y = 0). 


The stresses caused by IP were deduced 
in Section 2. For the two reflected waves that 
reach any point A situated in the interior of the 
rock mass, specific expressions are available. 

The reflected longitudinal wave (referred to by 
the subscript RP) has an angle of reflection equal 
to the angle of incidence, and its point of inci- 
dence is defined by the coordinates (x,, 0) where 
x, is easily found to be 


YaXe FYEXA 
YE TIA 


——— 


X | 


354 


and by the angle (see Fig. 5) 


Xp — X 

E I 
0, = arc tg — 
YE 


Now, the polar coordinates of A are (r,, 9, R 
$,) in relation to the point O, situated at the 
intersection of the normal to the free face pas- 
sing through E with the projection of the reflec- 
ted longitudinal wave. 

The stresses resulting from RP are then cal- 
culated with relation to O,, using the relation- 
ships from the theory of elasticity : (*) 


Ui (ri, 24) 


FE Um, (Tp) — 


rep = 0424) | 


| 


Vir To) à gã U,(rio 94, Lro) 
Cp E 
EU a (rs 9,) 
too = up e | SA | Uw, (Tip) — 
MA MRa 214 Top] 42346) 0 LHE 2seTro) 
Cp . 
JB. (r,,8,) 
E. pi 8 
Frgrp “x ll H o, Uw, (Tp) 
Top = "rep =—0 


The functions are calculated for the instant 


a, qq 
o 
Cp 
where r, = AO ,, and they have the expressions 


Dj (rio o Tr) = Ag (8,).AÇE Ui (T,o) 


rp 


Vilri, 8 To )=A,(9,).Bor ".Vwy( To) 
where A ,, is a function of /, that indicates the 
ratio between the amplitudes of RP and IP waves 
at the free face: 


A (0) sin 2 à, sin 2j masi Co/C yº cos” 2) 
pt egnçeos 
sin 2 ), sin 2) +( Cp/Cs y cos” 2) 
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and 


“e 
— sin 0, 
S 


sin j = — 


The wave shapes Uw, and Vw, are identical 


to the same functions of the incident wave, 
because the theory of reflection of elastic waves 


states the invariance of waves shapes with 
respect to reflections. 

Finally, let us consider the third wave encounte- 
ring the point A, which stresses are characterized 
by the subscript RS. The laws of reflection give 
for this case the position of the point of inci- 
dence defined by (0,x,) where 


X9<X, 


because C;< Cp. 
The conditions for the determination of x, as 
well as the angle 9, are: 


tg 3, = eai 


À “eadi 
sina, = — sin 6 
9 C 9 
S 


where 2, is the angle of incidence by means of 
which one stress wave departing from E is 
reflected in the direction of A, as a shear wave. 

Those three relationships constitue a system 
of 3 equations with 3 unknowns (2,, 2, and x,) 
that is solved in the computer program by means 
of a simple bracketing method. 

This system of polar coordinates will now 
be centered at O, (see Fig. 5) and each point in 
the plane has three components r,, 8,, 4. The 
time for arrival of the waves is 


and the corresponding stresses are given by: 


Ea Uol Tr) k U op (79,29) 


ro 08, 
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+ Uso (r,,0,) cotg O ] 


UmtT dU 
bw cá! K +92 3) dpfr 9189) do 


S 
bre 
rs d R 


+ LU op (179,45) cotg 1, 


ra Usw (Teo) k Usp (ra, 25) 
8, 


+(1+24) Usp (r,, 9%) cotg | 


To 


ú dU, (r , 
*rors = | op fes ta) Usw (T;) — 
dra 


Vo (r,, 0 ) 

ae Ti Vow (Te) — 

Usp (ro, 9,) 
E2 


Vw (TS | 


= 0 


"qtrs = Trtrs 


The displacement and particle velocity func- 
tions are in this case: 


h 
é SE 
Us (ro, GT) =A,d( IA ori E U,w (Tçs) 


-— 


rak 

Vo (1904 Tr) = Ap (0,)Botg "(=) Dari enad 
2 

where 


“2 sin 2i cos 29, 
A, (9) = É pq 


Eai 
sin 2i sin 29, + Eu cos” 29, 
Cs : 


sin i = — sin), 


and the following condition is valid 


U = U) = U,w 


w wl 
Vo ama Voo cm Vow 


meaning that the wave shape does not change on 
reflection. 


lhe next step in the stress calculation is to 
obtain the values to a common origin E. This 
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